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Self-assembly of polypeptides into fibrillar structures can be initiated by planar surfaces that 
interact favorably with certain residues. Using a coarse grained model, we systematically studied 
the folding and adsorption behavior of a /3-roll forming polypeptide. We find that there are two 
different folding pathways depending on the temperature: (i) at low temperature, the polypeptide 
folds in solution into a /3-roll before adsorbing onto the attractive surface, (ii) at higher tempera¬ 
ture, the polypeptide first adsorbs in a disordered state, and folds while on the surface. The folding 
temperature increases with increasing attraction, as the folded /3-roll is stabilized by the surface. 

Surprisingly, further increasing the attraction lowers the folding temperature again, as strong attrac¬ 
tion also stabilizes the adsorbed disordered state, which competes with folding of the polypeptide. 

Our results suggest that to enhance the folding, one should use a weakly attractive surface. They 
also explain the recent experimental observation of the non-monotonic effect of charge on the fibril 
formation on an oppositely charged surface [C. Charbonneau, et al, ACS Nano, 8, 2328 (2014)]. 

PACS numbers: 87.14.em, 87.15.A-, 87.15.lip, 87.15.Zg 


I. INTRODUCTION 


The interest in spontaneous fibril formation by struc¬ 
tural proteins derives not only from the link with neuro- 
degenerative diseases such as Alzheimers and Parkinsons, 
which have been traced to undesired amyloid fibril for¬ 
mation [T], but also from the fact that many natural 
and synthetic proteins form fibrils and hold promise for 
application as novel biomaterials HE]- In particular, the 
stimuli-responsive properties of fibrils have generated a 
strong interest in biomedical application [7]. Fibril for¬ 
mation is usually kinetically controlled and occurs via a 
nucleation-growth mechanism. While this mechanism is 
often believed to involve homogeneous nucleation, fibril 
formation can also be surface-induced. Such a mecha¬ 
nism has practical applications in industrial biosensors, 
in biotechnology, and in nanotechnology [3l |8]. Surface- 
mediated fibril formation occurs in oligopeptides (typi¬ 
cally around 20 residues long) [SHE] and amyloid beta 
peptide USHIS], but also for silk-elastine-like polymers 
nzmi]. The surface-induced fibril self-assembly process 
can occur via several routes. One involves a nucleation- 
and-growth process in solution, after which the pre¬ 
formed seed adsorbs on the surface, and continues to 
grow leHia ng n]. Another is by direct adsorption of 
the single molecules on the surface [9l [151 HZl fill- The 
morphology of surface-induced fibers depends on the pro¬ 
tein concentration, the physicochemical surface proper¬ 
ties, and environmental conditions such as temperature, 
pH, and ionic strength [3 [3 [TTl [El flBHIS] . 


We focus here on an example of designed biosynthetic 
peptide polymers based on silk-like and collagen-like 
sequences [5] • These polypeptides consist of three con¬ 
nected blocks, with the central one, inspired by sequences 
occurring in natural silk, a repeated octapeptide GAGA- 
GAGX, where G and A are glycine and alanine, respec¬ 
tively, and X is a polar residue such as glutamate or his¬ 
tidine. This middle block carries at either end a proline- 
rich and rather hydrophilic sequence inspired by natural 
collagen, which does not form any secondary structures 
but stays a random coil in aqueous solution. The silk-like 
block can fold into a /d-roll structure as soon as the charge 
on the polar residue is removed [H uni. Such /3-rolls 
have hydrophobic alanine rich faces, by which they self- 
assemble into long and filamentous stacks [Ellin]- At suf¬ 
ficiently high concentrations, these filaments form dilute 
hydrogels which are promising candidates as a matrix for 
artificial tissue [3 [6] • The question that concerns us here 
is how this process starts. Kinetic experiments indicated 
that individual molecules in solution do not readily fold 
under the conditions of the experiments [6l ITS] . and that 
fibril formation is governed by a nucleation-and-growth 
mechanism. In some cases (e.g., for X = histidine) ho¬ 
mogeneous nucleation seems to occur, but in others (e.g., 
when X = glutamate) this does not happen [5]. Recent 
experiments have shown that, in line with this, the pres¬ 
ence of a surface which is weakly attractive to certain 
residues in the peptide sequence can promote the for¬ 
mation of fibril structures, but the effect is rather subtle 
and the underlying physics remains unclear m- We have 
shown earlier that the formation of these fibrils is mainly 
triggered by the presence of a folded polypeptide, which 
then serves as a seed for the further growth of the long 
fibril [21]. The purpose of this work is to elucidate, by 
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using computer simulations, whether or not a flat surface 
can indeed take the role of the seed, and how the interac¬ 
tions between surface and polypeptide residues influence 
this. Since the experiments [18] suggested that the polar 
residue X is attracted towards the surface, we pay specific 
attention to this residue. To have an explicit model, we 
here choose X to be glutamate (E), but our conclusions 
should mutatis mutandis be valid for other choices. 


II. MODEL 

Although in principle atomistic models can provide in¬ 
sights into the folding of polypeptides, such all-atom sim¬ 
ulations are prohibitively expensive for the polymers con¬ 
sidered here. Therefore coarse grained modelling is the 
method of choice for making progress in understanding 
the physics of protein folding [20] [22H27] . Here, we em¬ 
ploy a coarse-grained polypeptide model in which each 
residue occupies a single site on a 3D cubic lattice, with 
all other sites considered as solvent Euim. In contrast 
to implicit solvent protein models, this model has been 
shown able to prevent artificial aggregation of proteins 
in their native state. Each residue has a unit vector in¬ 
dicating the direction of its side chain. Two essential 
features for describing folding, i.e., the formation of hy¬ 
drogen bonds and the directionality of side chains, are 
correctly captured in this highly efficient model. The 
total potential energy of the system is given hy E = 
^aa T ^solvent T ^hb T ^sterici where Eq^q^ and EgQiyQjii 
are are interaction potentials between residues, and be¬ 
tween a residue and solvent, respectively; the values of 
these have been obtained by comparison with experimen¬ 
tal data as shown in Table (all interaction potentials 
in this work are in reduced units). Ehij is the potential 
energy of formed hydrogen bonds, and Egteric represents 
the steric hindrance between consecutive residues in a 
polypeptide chain [29|. Two amino acids in contact in¬ 
teract only when their side chains are either in parallel 
or pointing towards each other. Similarly, interaction 
between a residue and solvent only exists when the side 
chain points to a solvent site. When a residue is not part 
of a turn in the backbone, it can adopt a strand confor¬ 
mation depending on the side chains. Two residues in 
contact and both in the strand state can form a hydro¬ 
gen bond with an energy Chb = —0.5, when their side 
chains are aligned. An energy penalty of = 0.55 is em¬ 
ployed to prevent the side chains of consecutive residues 
pointing in the same direction, thereby mimicking steric 
hindrance and restrictions in bond rotation. Due to its 
small size, alanine in a ^d-strand environment behaves as 
a considerably more hydrophobic residue than it actu¬ 
ally is; this effect is not well captured by the original 
parametrization of the potential. To compensate for this 
shortcoming, we vary the alanine-water interaction pa¬ 
rameter to investigate its influence. Configuration space 
is explored using a lattice Monte Carlo scheme with a 
classical set of moves m- 


TABLE I: Interaction matrix (in reduced units) for the 
residues in the used sequence [28]. Amino acids are denoted 
by their one-letter code (I = isoleucine, A = alanine, R = 
arginine, E = glutamate), and w and wall denote the solvent 
(water) and the residues on the surface. The hydrophobicity 
of alanine is varied via the A-w interaction {ca.w) as indicated. 

I A R E w wall 

I -0.79 -0.40 0.5 0.69 07 0 

A -0.34 0.49 0.77 [0.01, 0.6] 0 

R 0.43 -0.6 -0.57 0 

E 1.02 -0.78 [0.0, -3.0] 

w 0.0 0 

wall 0 



T 


FIG. 1: (color online) Heat capacity Cv — (E^) — (E)^ as a 
function of reduced temperature T* for various strength of 
attraction between the wall and glutamate CE.waii around the 
folding temperature for a single polypeptide with alanine hy¬ 
drophobicity eA,w = 0.6. Inset: (Top) Heat capacity Cy as 
a function of T* for the whole temperature range including 
both folding and adsorption of the polypeptide. (Bottom) 
The free energy barrier height of peptide folding, AF*/ksT, 
as a function of surface attraction strength \eE,waii\ at a con¬ 
stant temperature T* = 0.43. 


As lattice models in general are not expected to fold 
natural sequences, we designed a sequence to fold into 
the desired /3-roll structure. As explained in Ref[2T] 
proper folding of a /3-roll on a cubic lattice demands 
a palindromic sequence and a anti-parallel sheet topol¬ 
ogy. We restrict the silk part to 80 residues to make 
the calculation tractable, for which the design procedure 
yielded the optimal sequence {E{AI)sRE{IA)sR)6- The 
replacement of glycine with isoleucine is not unrealistic, 
as glycines in ^d-sheets are more hydrophobic than the av¬ 
erage glycine [21] . The extra (arginine) residue has been 
introduced to fit the /3-roll structure on the lattice. 

We performed replica exchange Monte Carlo (REMC) 
simulations for one polypeptide chain with sequence 
{E{AI)^RE{IA)^R)q in a simulation box of 100 x 100 x 
100 lattice sites with periodic boundary conditions in x 
and y directions. To mimic the experimental situation, 
in which the X (here, E) residue is electrostatically at- 
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FIG. 2: (color online) Top: Free energy landscapes as a function of the distance between the polypeptide center of mass and the 
wall, Rwc, and the total number of hydrogen bonds, Htot, for a single polypeptide with the alanine hydrophobicity eA,w = 0.6, 
and wall attraction es^waii = —0.6 at different temperatures. Insets: Snapshots of a typical configuration corresponding to the 
local minima on the free energy landscapes indicated by the arrows (from left to right): an unfolded molecule, a folded molecule 
at the surface, and a folded one in solution. Residue color coding: alanine(A)=yellow, isoleucine(I)=white, glutamate(E)=red, 
arginine(R)=blue, and wall=green. Bottom: Free energy landscape as a function of Rwc and Htot for the system containing a 
single polypeptide with eA,w = 0.6, and various wall attraction eE,waii at the temperature T* = 0.451. 


tracted to the surface we put two parallel walls on op¬ 
posite sides of the box in the z direction consisting of 
sites which exclusively attract E with a variable strength 
^E,waii (see Table [^. We ensure that the size of the sim¬ 
ulation box is large enough not to influence the folding 
of the polypeptide. Each REMC simulation consisted 
of 48 replicas with a (reduced) temperature distribution 
around the transition temperature, which is optimized by 
a feedback-optimization algorithm m- For each replica, 
we performed 6 x 10^^ MC cycles with an exchange at¬ 
tempt every 1000 cycles. The first 10^^ moves were for 
equilibration. Employing the virtual-move parallel tem¬ 
pering method during the production m ensured opti¬ 
mal use of the simulation data. 


III. RESULTS 

We first simulate the folding of the polypeptide with 
a relatively strong alanine hydrophobicity eA,w = 0.6, 
which is close to the experimental situation [21]. The 
heat capacity Cy calculated from the energy fluctua¬ 
tions, as a function of the reduced temperature T* is 
shown in Fig. for different surface attraction strengths. 
Even when the surface is not attractive at all, i.e. at 
^E,waii = 0, the folding temperature is somewhat higher 


than the bulk value, implying that the surface stabi¬ 
lizes the folded structure. This stability follows from 
the hydrophobicity of the outside surface of the folded 
/3-roll structure (exposing mostly alanine), which prefers 
the surface. With increasing attraction strength, from 
\^E,waii\ = 0 to 2.0, the folding temperature rises. The 
corresponding folding free energy barrier height is shown 
in the inset of Fig. for a temperature T* = 0.43 at 
which the polypeptide cannot spontaneously fold in the 
bulk solution. The barrier dramatically decreases with 
attraction to the wall, indicating the attraction promotes 
the folding of the polypeptide. 

Fig-i top shows typical free energy landscapes for 
^E,waii = —0.6 at different temperatures, as a function 
of the distance between the mass center of the polypep¬ 
tide and the surface, Rwc^ and the total number of hy¬ 
drogen bonds formed in the system, Htot- At low tem¬ 
perature T* = 0.35 the stable phase of the system is 
the folded /3-roll structure, represented by a global min¬ 
imum at {Rwc ~ ~ 30). Another local mini¬ 

mum at {Rwc ~ 30, Htot ~ 20) corresponds to a folded 
/3-roll floating at a distance from the surface. At this 
temperature the polypeptide is likely to first fold in so¬ 
lution before adsorbing onto the attractive surface. At 
T* = 0.42961 (Fig. top left), just above the folding 
transition, the global minimum is at {Rwc ~ 3, Htot ~ 0), 
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corresponding to a disordered polypeptide adsorbed on 
the surface. A channel appears that connects the dis¬ 
ordered and the folded structures on the surface. At 
T* = 0.41115, slightly below the wall-induced folding 
temperature Tf = 0.4166, but still above the bulk fold¬ 
ing temperature = 0.406 (Fig.j^top middle panel), 
the adsorbed /3-roll structure becomes the stable state 
of the system. This strongly suggests that the adsorbed 
disordered polypeptide folds while in contact with the 
surface. The presence of the attractive surface thus not 
only stabilizes the /3-roll, but, as a first step towards fold¬ 
ing, also brings the disordered structure towards the sur¬ 
face. For this particular /3-roll structure, only half of the 
octapeptide strands, and hence half of the E residues, 
are in direct contact with the surface. This is similar to 
the structure found in atomistic simulations where the E 
residues were pointing in opposite directions on the two 
sides of the folded ;d-roll [2^ . 

Upon increasing the surface attraction strength 
\^E,waii\ = 2.0 to 3.0, intriguingly, the folding temper¬ 
ature moves down, and the free energy barrier for folding 
increases again. Simultaneously, a new high temperature 
peak appears in the heat capacity. This peak corresponds 
to the adsorption/desorption transition of the polypep¬ 
tide, and shifts to higher temperatures with increasing 
surface attraction strength due to the enhanced stability 
of the adsorbed phase. To see why the folding temper¬ 
ature moves down, we again plot free energy landscapes 
as a function of R^c and Htot ia Fig. (bottom) for a 
fixed temperature (T* = 0.451) above the folding tran¬ 
sition for varying surface attractions from \eE,waii\ = 2.0 
to 3.0. In these diagrams the two local minima, located 
at (i?^c ~ l.b.Htot ~ 20) and {R^c - - 0), 

correspond to the adsorbed folded and disordered state, 
respectively. The dissolved, non-adsorbed state is now 
completely suppressed. While at eE,waii = —2.0 the two 
states differ only slightly in free energy, increasing the at¬ 
traction stabilizes the disordered state, and enlarges this 
difference to as much as 10k bT for eE,waii = —3.0. Along 
with this change, the average peptide-wall distance R^c 
decreases as the chain flattens out on the surface. As 
in the folded structure only one half of the glutamate 
residues are in direct contact with the surface, the re¬ 
maining glutamate residues also tend to stick to the sur¬ 
face when the attraction becomes too strong, but this can 
only happen if the /3-roll structure is destroyed. Such be¬ 
havior is also observed in simulations of small proteins 
next to an attractive wall [321436] . Hence, the effect of 
attraction on the efficiency of surface-induced folding is 
non-monotonic, because of the specific topology of the 
folded structure. We note that as the surface adsorbed 
P-ToW exposes one alanine-rich face toward the solution, 
it likely will act as a nucleus for fibril growth by favorable 
interaction between two such faces [21] . 

The folding mechanism of the polypeptide on an at¬ 
tractive surface is explained with schematic plots of the 
free energy vs. temperature in Eig. [^-d. Eor low sur¬ 
face attraction, the polypeptide folds in the bulk solution 




T T 



FIG. 3: (color online) (a-d) The schematic illustration of the 
free energies of three states of the system F as a function 
of temperature T with increasing wall attraction from (a) to 
(d). Black, red, and blue lines are the free energies for folded 
state, disorder state on the wall and the disorder state in the 
bulk solution, respectively, and Tf and Ta are the folding and 
adsorption temperatures, respectively, (e) Critical tempera¬ 
ture, Tc, for the transitions, i.e. folding (filled symbols) and 
adsorption (open symbols), as a function of the wall attraction 
\eE,waii\ for the polypeptide with various alanine hydropho- 
bicities eA,w The horizontal dashed lines indicate the folding 
temperature of polypeptide in bulk solution. 


when decreasing the temperature of the system, which 
occurs as a direct transition from the disordered state 
to folded state at Tj (Fig. |^). With increasing surface 
attraction, the folded state is stabilized, the free energy 
of the folded state decreases, and the folding tempera¬ 
ture decreases (Fig. If)- A very strong surface attraction 
also stabilizes a disordered polypeptide on the surface, 
which splits the folding of polypeptide into two steps: 
with decreasing temperature, the disordered polypeptide 
first adsorbs onto the surface at temperature T^, before 
folding on the surface at Tj (Fig. [^). Stronger surface 
attraction increases the adsorption temperature Ta and 
lowers the folding temperature Tj (Fig.j^). Fig.plots 
the two transition temperatures Ta and Tj, as a func¬ 
tion of surface attraction at various alanine hydropho- 
bicities For all three cases, the surface effect is 

non-monotonic, with Tf showing a maximum around 
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^E,waii = “2 (lower curves) while the adsorption tem¬ 
peratures increase monotonically (upper curves). The 
mechanism seems to be independent of the hydrophobic- 
ity of the outer face, which is crucial for the stacking of 
the /3-roll [21]. We note that this explanation implies 
that there is a certain value for eA,w at which the dis¬ 
ordered polypeptide just barely becomes stable. This is 
similar to the situation in which a liquid phase becomes 
stable, as a function of the pressure in a simple liquid 
phase diagram. While in principle this special value for 
eA,w does exist, it is hard to measure, even in our lat¬ 
tice model, due to the large error bars connected to the 
estimating the transition temperature. 

IV. CONCLUSIONS 

In conclusion, by performing computer simulations of 
a coarse-grained lattice model polypeptide, we obtain 
mechanistic insights into the effect of surface attraction 
on the folding of a single disordered polypeptide with a 
silk-like sequence. We find that increasing attraction be¬ 
tween surface and polar residue stabilizes folding, i.e. the 
folding temperature increases, and the free energy bar¬ 
rier for the folding of polypeptide decreases. A pathway 
for surface-induced folding opens up along which the sur¬ 


face first captures the disordered chain which then can 
pass over to the folded state, even above the bulk folding 
temperature. In contrast, at lower temperature, folding 
can first occur in solution after which the folded struc¬ 
ture adsorbs to the wall. At strong surface attraction 
the polypeptide tends to flatten out entirely, rendering 
the folded structure unstable with respect to the disor¬ 
dered state, with a higher folding barrier. 

Our results suggest, therefore, that in order to pro¬ 
mote surface-induced folding to seed the hierarchical self- 
assembly of protein fibrils, one should experimentally op¬ 
erate in a window of relatively weak binding where the 
chain flattening does not occur. We stress that these 
findings are generic, and should translate to other protein 
systems. 
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